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Abstract 

We review realistic models that reproduce quantum theory in some limit and yield poten- 
tially new physics outside that limit. In particular, we consider deterministic hidden-variables 
theories (such as the pilot-wave model) and their extension to "quantum nonequilibrium," and 
we consider the continuous spontaneous localization model of wave function collapse. Other 
models are briefly discussed. 
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1 Introduction 



According to the so-called "Copenhagen Interpretation," standard quantum theory is limited 

to describing experimental sitnations. It is at once remarkably successful in its predictions, and 
remarkably ill-defined in its conceptual structure: what is an experiment? what physical objects 
do or do not require quantization? how are the states realized in nature to be characterized? 
how and when is the wave-function "collapse postulate" to be invoked? Because of its siicccss, 
one may suspect that quantum theory can be promoted from a theory of measurement to a 
theory of reality. But, that requires there to be an unambiguous specification (S) of the possible 
real states of Nature and their probabilities of being realized. 

There are several approaches that attempt to achieve S. The more conservative approaches 
(consistent histories, environmental decoherence, many worlds) do not produce any predictions 
that differ from the standard ones because they do not tamper with the usual basic mathemat- 
ical formalism. Rather, they utilize structures compatible with standard quantum theory to 
elucidate S. These approaches, which will not be discussed in this article, have arguably been 
less successful so far at achieving S than approaches that introduce significant alterations to 
quantum theory. 

This article will largely deal with the two most well-developed realistic models that reproduce 

quantum theory in some limit and yield potentially new and testable physics outside that limit. 
First, the pilot-wave model, which will be discussed in the broader context of "hidden-variables 
theories." Second, the continuous spontaneous localization (CSL) model, which describes wave 
function collapse as a physical process. Other related models will also be discussed briefly. 

Due to bibliographic space limitations, this article contains a number of uncited references, 
of the form "(author) in (year)." Those in the next section can be found in Valentini (2002b, 
2004a, b) or at www.arxiv.org. Those in subsequent sections can be found in Adler (2004), Bassi 
and Ghirardi (2003), Pearle (1999) (or in subsequent papers by these authors, or directly, at 
www.arxiv.org), and in Wallstrom (1994). 

2 Hidden Variables and Quantum Nonequilibrium 

A deterministic hidden- variables theory defines a mapping ui = uj{M, A) from initial hidden 

parameters A (defined, e.g., at the time of preparation of a quantum state) to final outcomes 
u of quantum measurements. The mapping depends on macroscopic experimental settings M, 
and fixes the outcome for each run of the experiment. Bell's theorem of 1964 shows that, for 
entangled quantum states of widely separated systems, the mapping must be nonlocal: some 
outcomes for (at least) one system must depend on the setting for another distant system. 

In a viable theory, the statistics of quantum measurement outcomes - over an ensemble of 
experimental trials with fixed settings M will agree with quantum theory for some special 
distribution pqt(A) of hidden variables. For example, expectation values will coincide with the 
predictions of the Born rule 

(a;)QT = J dX pqt(A)w(M,A) = Tr(pfi) 

for an appropriate density operator p and Hermitian observable Q. (As is customary in this 
context, J dX is to be understood as a generalised sum.) However, given the mapping lo = 
u){M, A) for individual trials, one may, in principle, consider nonstandard distributions p(A) ^ 
Pqt(A) that yield statistics outside the domain of ordinary quantum theory (Valentini 1991, 
2002a) . We may say that such distributions correspond to a state of quantum nonequilibrium. 

Quantum nonequilibrium is characterised by the breakdown of a number of basic quantum 
constraints. In particular, nonlocal signals appear at the statistical level. We shall first illustrate 
this for the hidden-variables model of de Broglie and Bohm. Then we shall generalize the 
discussion to all (deterministic) hidden-variables theories. 
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At present there is no experimental evidence for quantum nonequilibrium in nature. How- 
ever, from a hidden-variables perspective, it is natural to explore the theoretical properties of 
nonequilibrium distributions, and to search experimentally for the statistical anomalies associ- 
ated with them. 

From this point of view, quantum theory is a special case of a wider physics, much as thermal 
physics is a special case of a wider (nonequilibrium) physics. (The special distribution /3qt(A) 
is analogous to, say. Maxwell's distribution of molecular speeds.) Quantum physics may be 
compared with the physics of global thermal equilibrium, which is characterised by constraints - 
such as the impossibility of converting heat into work (in the absence of temperature differences) 
- that are not fundamental but contingent on the state. Similarly, quantum constraints such as 
statistical locality (the impossibility of converting entanglement into a practical signal) are seen 
as contingencies of pqt{\). 



2.1 Pilot- Wave Theory 

The de Broglie-Bohm "pilot-wave theory" - as it was originally called by de Broglie, who first 
presented it at the Fifth Solvay Congress in 1927 - is the classic example of a deterministic 
hidden- variables theory of broad scope (Bohm 1952; Bell 1987; Holland 1993). We shall use it 
to illustrate the above ideas. Later, the discussion will be generalised to arbitrary theories. 

In pilot- wave dynamics, an individual closed system with (configuration-space) wave function 
^{X,t) satisfying the Schrodinger equation 

a* 

^n-^H^ (1) 
has an actual configuration X{t) with velocity 

Mt) - (2) 



where J — J [^j — J{X, t) satisfies the continuity equation 

dt = ' 

(which follows from Q). In quantum theory, J is the "probability current" . In pilot- wave theory, 
5* is an objective physical field (on configuration space) guiding the motion of an individual 
system. 

Here, the objective state (or ontology) for a closed system is given by 5" and X. A probability 
distribution for X - discussed below - completes an unambiguous specification S (as mentioned 
in the introduction). 

Pilot-wave dynamics may be applied to any quantum system with a locally conserved current 
in configuration space. Thus, X may represent a many-body system, or the configuration of a 
continuous field, or perhaps some other entity. 

For example, at low energies, for a system of N particles with positions Xi(i) and masses rrij 
[i — 1, 2, N), with an external potential V , Q (with X = (xi, X2, xjv)) reads 

^h^^y^-^^l^ + V^ (4) 
dt ^ 2m., ' ^ ' 



while (O has components 
(where = \^\e^^l^^^). 



^ = f Zi^) . M (5) 

dt TO,- V * / TO,- 
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In general, ^ and ^ determine X{t) for an individual system, given the initial condtions 
X{Q), ^{X, 0) at t = 0. For an arbitrary initial distribution P{X, 0), over an ensemble with the 
same wave function \E'(X, 0), the evolution P{X, t) of the distribution is given by the continuity 
equation 

dP 

— + V • {PX) - (6) 

The outcome of an experiment is determined by X{Q), ^'(X, 0), which may be identified with 
A. For an ensemble with the same ^'(X, 0), we have A = -'^"(O). 

2.1.1 Quantum Equilibrium 

From © and ©, if we assume P(X,0) = 0)|^ at i = 0, we obtain P[X,t) = - 

the Born-rule distribution of configurations - at all times t. 

Quantum measurements are, like any other process, described and explained in terms of 
evolving configurations. For measurement devices whose pointer readings reduce to configu- 
rations, the distribution of outcomes of quantum measurements will match the statistical pre- 
dictions of quantum theory (Bohm 1952; Bell 1987; Diirr, Goldstein and Zanghi 2003). Thus 
quantum theory emerges phenomcnologically for a "quantum equilibrium" ensemble with dis- 
tribution P{X,t) = (or p{\) = pqt(A)). 

2.1.2 Quantum Nonequilibrium 

In principle, as we saw for general hidden-variables theories, we may consider a nonequilibrium 
distribution P{X^ 0) ^ |^(^, 0)|^ of initial configurations while retaining the same deterministic 
dynamics (0, 10) for individual systems (Valentini 1991). The time evolution of P{X,t) will be 
determined by 

As we shall see, in appropriate circumstances (with a sufficiently complicated velocity field 
X), © generates relaxation P ^ jvPI on a coarse-grained level, much as the analogous classical 
evolution on phase space generates thermal relaxation. But for as long as the ensemble is in 
nonequilibrium, the statistics of outcomes of quantum measurements will disagree with quantum 
theory. 

Quantum nonequilibrium may have existed in the very early universe, with relaxation to 
equilibrium occurring soon after the big bang. Thus, a hidden-variables analogue of the clas- 
sical thermodynamic "heat death of the universe" may have actually taken place (Valentini 
1991). Even so, relic cosmological particles that decoupled sufficiently early could still be in 
nonequilibrium today, as suggested by Valentini in 1996 and 2001. It has also been speculated 
that nonequilibrium could be generated in systems entangled with degrees of freedom behind a 
black- hole event horizon (Valentini 2004a). 

Experimental searches for nonequilibrium have been proposed. Nonequilibrium could be 
detected by the statistical analysis of random samples of particles taken from a parent population 
of (for example) relics from the early universe. Once the parent distribution is known, the rest 
of the population could be used as a resource, to perform tasks that are currently impossible 
(Valentini 2002b). 

2.2 iif-Theorem: Relaxation to Equilibrium 

Before discussing the potential uses of nonequilibrium, we should first explain why all systems 
probed so far have been found in the equilibrium state P = l^-l . This distribution may be 
accounted for along the lines of classical statistical mechanics, noting that all currently-accessible 
systems have had a long and violent astrophysical history. 

Dividing configuration space into small cells, and introducing coarse-grained quantities P, 
l^*!^, a general argument for relaxation P ^ \'^\^ is based on an analog of the classical coarse- 
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graining H-theoiem. The coarse-grained ff-function 



(7) 



(minus the relative entropy of P with respect to l^*!^) obeys the /^-theorem (Valentini 1991) 

H{t) < H{0) 

(assuming no initial fine-grained microstructure in P and l^*]^). Here, H > for all P, I'^f and 
H = ii and only if P = everywhere. 

The if-theorem expresses the fact that P and | \1/ 1 behave like two "fluids" that are "stirred" 
by the same velocity field X, so that P and \^\ tend to become indistinguishable on a coarse- 
grained level. Like its classical analog, the theorem provides a general understanding of how 
equilibrium is approached, while not proving that equilibrium is actually reached. (And of 
course, for some simple systems - such as a particle in the ground state of a box, for which 
the velocity field VS/m. vanishes there is no relaxation at all.) A strict decrease of H{t) 
immediately after f = is guaranteed if Xq ■ V(Po/ 1^*0 1 ) has non-zero spatial variance over a 
coarse-graining cell, as shown by Valentini in 1992 and 2001. 

A relaxation timescalc t may be defined by l/r^ = - H / dt:^) Hq. For a single par- 
ticle with quantum energy spread AE, a crude estimate given by Valentini in 2001 yields 
T ~ (l/£);iVw^/2(AP)3/2, where e is the coarse-graining length. For wave functions that 
are superpositions of many energy eigenfunctions, the velocity field (generally) varic;s rapidly, 
and detailed numerical simulations (in two dimensions) show that relaxation occurs with an 
approximately exponential decay H{t) w .ffoe~*^*°, with a time constant tc of order r (Valentini 
and Westman 2005). 

Equilibrium is then to be expected for particles emerging from the violence of the big bang. 
The possibility is still open that relics from very early times may not have reached equilibrium 
before decoupling. 



2.3 Nonlocal Signaling 

We now show how nonequilibrium, if it were ever discovered, could be used for nonlocal signaling. 

Pilot-wave dynamics is nonlocal. For a pair of particles A, B with entangled wave function 
4'(x^,XB,t), the velocity 'k.A{t) = VyiS'(Xyii, xs, t)/TOyi of A depends instantaneously on x^, 
and local operations at B - such as switching on a potential - instantaneously affect the motion 
of A. For an ensemble P(x^,xs,t) = |^(xA,xs,f)p, local operations at B have no statistical 
effect at A: the individual nonlocal effects vanish upon averaging over an equilibrium ensemble. 

Nonlocality is (generally) hidden by statistical noise only in quantum equilibrium. If instead 
P(x^,xb,0) 7^ |'J'(x^,xb,0)|2, a local change in the Hamiltonian at B generally induces an 
instantaneous change in the marginal pyi(xyi, f) = J d'^xs P(x/i,XB,t) at A. For example, in 
one dimension a sudden change Hb H'g in the Hamiltonian at B induces a change ApA = 
PA{xA,t) — pa{xa,0) (for small t) (Valentini 1991), 

A ^ f i \ f ^ ^ P{xa,xb,Q)-\^{xa,xbM'' \ ic,. 

= -^d^A r^'U '^^'^ m^A,XB,0)\^ ) 

(Here niA = tub = m, a{xA) depends on ^'(a;^, x^, O), while b{xB) also depends on H'^ and 
vanishes if H'g = Hb-) The signal is generally nonzero if Pq 7^ l^'ol^- 

Nonlocal signals do not lead to causal paradoxes if, at the hidden-variable level, there is 
a preferred foliation of spacetime with a time parameter that defines a fundamental causal 
sequence. Such signals, if they were observed, would define an absolute simultaneity as discussed 
by Valentini in 1992 and 2005. Note that in pilot- wave field theory, Lorentz invariance emerges as 
a phenomenological symmetry of the equilibrium state, conditional on the structure of the field- 
theoretical Hamiltonian (as discussed by Bohm and Hiley in 1984, Bohm, Hiley and Kaloyerou 
in 1987, and Valentini in 1992 and 1996). 
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2.4 Subquantum Measurement 

In principle, nonequilibrium particles could also be used to perform "subquantum measure- 
ments" on ordinary, equilibrium systems. We illustrate this with an exactly solvable one- 
dimensional model (Valentini 2002b). 

Consider an apparatus "pointer" coordinate y. with known wave function go{y) and known 
(ensemble) distribution 7ro(j/) 7^ |<7o(y)l^i where 7ro(?/) has been deduced by statistical analysis 
of random samples from a parent population with known wave function go{y). (We assume that 
relaxation may be neglected: for example, if go is a box ground state, y = and 7ro(?/) is static.) 
Consider also a "system" coordinate x with known wave function "tpoix) and known distribution 
Po{x) = \il}o{x)\ . If 7ro(y) is arbitrarily narrow, xq can be measured without disturbing V'o(a^)) 
to arbitrary accuracy (violating the uncertainty principle). 

To do this, at t = we switch on an interaction Hamiltonian H = axpy, where a is a constant 
and Py is canonically conjugate to y. For relatively large a, we may neglect the Hamiltonians 
of x and y. For ^ = ^{x,y,t), we then have d'^/dt = —axd^/dy. For we have the 
continuity equation d \'^\^ /dt = —axd \'^\^ /dy, which implies the hidden-variable velocity fields 
a; = 0, y = ax and trajectories x{t) = Xq, y{t) = j/o + ax^t. 

The initial product ^o{x,y) = ipo{x)go{y) evolves into ^{x,y,t) = ipo{x)go{y — axt). For 
at ^ Q (with a large but fixed), ^i{x,y,t) ipo(x)go{y) and ipo{x) is imdisturbed: for small 
at, a standard quantum pointer with the coordinate y would yield negligible information about 
xq. Yet, for arbitrarily small at, the hidden- variable pointer coordinate y{t) = yo + axot does 
contain complete information about xq (and x(t) = xq). This "subquantum" information will 
be visible to us if no{y) is sufficiently narrow. 

For, over an ensemble of similar experiments, with initial joint distribution Po{x,y) = 
\il>o{x)\^ TTo{y) (equilibrium for x and nonequilibrium for y), the continuity equation dP/dt = 
—axdP/dy imphes that P{x,y,t) = |V'o(a;)| 7ro(j/ — axt). If 7ro(j/) is localised around y = 
("■0(2/) = for \y\ > w/2), then a standard (faithful) measurement of y with result j/moas will 
imply that x lies in the interval {ymcas/at — w/2at, ynicas/at + w/2at) (so that P{x, y, t) ^ 0). 
Taking the simultaneous limits at 0, w ^ 0, with w/at — » 0, the midpoint ymeas/at xq 
(since j/mcas = yo + CLXot and \yo\ < w/2), while the error w/2at 0. 

If w is arbitrarily small, a sequence of such measurements will determine the hidden trajec- 
tory x{t) without disturbing tp{x,t), to arbitrary accuracy. 

2.5 Subquantum Information and Computation 

From a hidden- variables perspective, immense physical resources are hidden from us by equilib- 
rium statistical noise. Quantum nonequilibrium would probably be as useful technologically as 
thermal or chemical nonequilibrium. 

2.5.1 Distinguishing Nonorthogonal States 

In quantum theory, nonorthogonal states \tpi), IV'2) ((V'ilV'2) 7^ 0) cannot be distinguished 
without disturbing them. This theorem breaks down in quantum nonequilibrium (Valentini 
2002b). For example, if \tpi), IV'2) are distinct states of a single spinless particle, then the 
associated de Broglie-Bohm velocity fields will in general be different, even if (V'ilV'2) 7^ 0, and 
so will the hidden- variable trajectories. Subquantum measurement of the trajectories could then 
distinguish the states \tpi), IV'2)- 

2.5.2 Breaking Quantum Cryptography 

The security of standard protocols for quantum key distribution depends on the validity of 
the laws of quantum theory. These protocols would become insecure given the availability of 
nonequilibrium systems (Valentini 2002b). 
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The protocols known as BB84 and B92 depend on the impossibility of distinguishing nonorthog- 
onal quantum states without disturbing them. An eavesdropper in possession of nonequilibrium 
particles could distinguish the nonorthogonal states being transmitted between two parties, and 
so read the supposedly secret key. Further, if subquantum measurements allow an eavesdropper 
to predict quantum measurement outcomes at each "wing" of a (bipartite) entangled state, then 
the EPR (Einstein-Podolsky-Rosen) protocol also becomes insecure. 

2.5.3 Subquantum Computation 

It has been suggested that nonequilibrium physics would be computationally more powerful than 
quantum theory, because of the ability to distinguish nonorthogonal states (Valentini 2002b). 
However, this ability depends on the (less-than-quantum) dispersion w of the nonequilibrium 
ensemble. A well-defined model of compiitational complexity requires that the resources be 
quantified in some way. Here, a key question is how the required w scales with the size of the 
computational task. So far, no rigorous results are known. 

2.6 Extension to All Deterministic Hidden- Vciriables Theories 

Let us now discuss arbitrary (deterministic) theories. 

2.6.1 Nonlocal Signaling 

Consider a pair of two-state quantum systems A and B, which are widely separated and in the 
singlet state. Quantum measurements of observables a a = • &a, = ms • (Ts (where 
m^, ms are unit vectors in Bloch space and a a, o's are Pauli spin operators) yield outcomes 
o'A, = ±1, in the ratio 1 : 1 at each wing, with a correlation (o'ao'b) = — • niB. Bell's 
theorem shows that for a hidden- variables theory to reproduce this correlation - upon averaging 
over an equilibrium ensemble with distribution pqt (A) - it must take the nonlocal form 

= o-A(m^,mB, A), (7s = CTB(m^,mB, A) (9) 

More precisely, to obtain {aAcrB)Qj^ = —^A ■ niB (where (ctactbIqt — / ^A PqtW'^ao'b), at 
least one of cr^, ctb must depend on the measurement setting at the distant wing. Without loss 
of generality, we assume that a a depends on niB. 

For an arbitrary nonequilibrium ensemble with distribution p(A) ^ pqt(A), in general 
(cr^ifTB) = J dX p{X)aAcrB difi'ers from — • mB, and the outcomes cr^, ctb = ±1 occur 
in a ratio different from 1:1. Further, a change of setting mB ^ at B will generally induce 
a change in the outcome statistics at A, yielding a nonlocal signal at the statistical level. To 
see this, note that, in a nonlocal theory, the "transition sets" 

Ta{-,+) = {A|crA(mA,mB, A) = -1, crA(m^, m'^. A) = +1} 
Ta{+,-) = {A|crA(mA,mB, A) = +1, aA{mA,m.'s, X) = -1} 

cannot be empty for arbitrary settings. Yet, in quantum equilibrium, the outcomes aA = ±1 
occur in the ratio 1 : 1 for all settings, so the transition sets must have equal equilibrium measure, 
Mqt[7a(— ,+)] = A^Qt[7a(+, — )] {dfJ-QT = PQT(A)dA). That is, the fraction of the equilibrium 
ensemble making the transition (7^ = — 1^(ta = +1 under hib ^ must equal the fraction 
making the reverse transition aA = +^ cta = —1- (This "detailed balancing" is analogous 
to the principle of detailed balance in statistical mechanics.) Since Ta(— ,+), Ta{+,—) are 
fixed by the deterministic mapping, they are independent of the ensemble distribution /5(A). 
Thus, for p(A) ^ Pqt(A), in general /x[Ta(— , +)] ^ ij[Ta{+, — )] {dp = p{X)dX): the fraction of 
the nonequilibrium ensemble making the transition (Ta = — 1— *cr^ = +l will not in general 
balance the fraction making the reverse transition. The outcome ratio at A will then change 
under mB ^ m^ and there will be an instantaneous signal at the statistical level from B to A 
(Valentini 2002a). 
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Thus, in any deterministic hidden-variables theory, nonequilibrium distributions p{X) ^ 
Pqt(A) generally allow entanglement to be used for nonlocal signaling (just as, in ordinary 
statistical physics, differences of temperature make it possible to convert heat into work). 

2.6.2 Experimental Signature of Nonequilibrium 



relates statistics from distinct, "incompatible" experiments. In nonequilibrium, such additivity 
generically breaks down (Valentini 2004b). 

Further, for a two-state system with observables m • &, the "dot-product" structure of the 
quantum expectation (m ■ a) =Tr (pm • it) = m • P (for some Bloch vector P) is equivalent to 
expectation additivity (Valentini 2004b). Nonadditive expectations then provide a convenient 
signature of nonequilibrium for any two-state system. For example, the sinusoidal modulation 
of the quantum transmission probability for a single photon through a polariser 



(where an angle 9 on the Bloch sphere corresponds to a physical angle 8 = 6*72) will generically 
break down in nonequilibrium. Deviations from H10() would provide an unambiguous violation 
of quantum theory (Valentini 2004b). 

Such deviations were searched for by Papaliolios in 1967, using laboratory photons and 
successive polarisation measurements over very short times, to test a hidden-variables theory 
(distinct from pilot- wave theory) due to Bohm and Bub (1966), in which quantum measure- 
ments generate nonequilibrium for short times. Experimentally, successive measurements over 
timescales ~ 10"^'^ sec agreed with the (quantum) sinusoidal modulation cos^ to < 1%. Sim- 
ilar tests might be performed with photons of a more exotic origin. 

3 Continuous Spontaneous Localization Model (CSL) 

The basic postulate of CSL is that the state vector \tp,t) represents reality. Since, for example, 
in describing a measurement, the usual Schrodinger evolution readily takes a real state into a 
non-real state, that is, into a superposition of real states (such as apparatus states describing 
different experimental outcomes), CSL requires a modification of Schrodinger's evolution. To 
the Hamiltonian is added a term which depends upon a classical randomly fluctuating field 
w(x, i) and a mass-density operator A(x, i). This term acts to collapse a superposition of 
states, which differ in their spatial distribution of mass density, to one of these states. The 
rate of collapse is very slow for a superposition involving a few particles, but very fast for a 
superposition of macroscopically different states. Thus, very rapidly, what you see (in nature) 
is what you get (from the theory). Each state vector evolving under each ^(x, t) corresponds to 
a realizable state, and a rule is given for how to associate a probability with each. In this way, 
an unambiguous specification S as mentioned in the introduction is achieved. 

3.1 Requirements for Stochastic Collapse Dynamics 

Consider a normalized state vector \ip,t) — X^n '^n(^)l'^n) (('^nkn') — ^rm') which undergoes a 
stochastic dynamical collapse process. This means that, starting from the initial superposition 
at t = 0, for each run of the process, the squared amplitudes Xn{t) = |a„(t)p fluctuate until all 
but one vanish, that is, x„i(oo) — 1, (a;^m(cxD) — 0) with probability Xm(0). 

This may be achieved simply, assuming negligible effect of the usual Schrodinger evolution. 




PQT(e) = - (1 + (m • ^)) = - (1 + Pcos2e) 



(10) 
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if the stochastic process enjoys the following properties (Pearle 1979): 



^x„(<) = l (11a) 

n 

'^^^XniO) (lib) 

Xn{oo)xm{oo) — for TO ^ n, (11c) 



where the overbar indicates the ensemble average at the indicated time. The only way that a sum 
of products of non-negative terms can vanish is for at least one term in each product to vanish. 
Thus, according to (|llcl) . for each run, at least one of each pair {a;„(cxD), 2;„j(oo)} (n ^ to) 
must vanish. This means that at most one Xn{oo) might not vanish and, by Hlla|l . applied 
at t = oo, one x„(oo) must not vanish and, in fact, must equal 1: hence, each run produces 
collapse. Now, let the probability of the outcome {x„(oo) = 1, a;^„(oo) = 0} be denoted P„. 
Since a;„(oo) = 1 • P„ + X^m^^n ' ^™ ~ -f" Iben, according to the Martingale property (|llb|l . 
applied a,t t = oo, Pn — Xn{0): hence the ensemble of runs produces the probability postulated 
by the usual "collapse rule" of standard quantum theory. 

A (nonquantum) stochastic process which obeys these equations is the gambler's ruin game. 
Suppose one gambler initially possesses the fraction a;i(0) of their joint wealth, and the other 
has the fraction 2:2(0). They toss a coin: heads, a dollar goes from gambler 1 to gambler 2, tails 
the dollar goes the other way. (|lla|l is satisfied since the sum of money in the game remains 
constant, Ijllbl) holds because it is a fair game, and Ullcfl holds because each game eventually 
ends. Thus, gambler i wins all the money with probability Xi(0). 

3.2 CSL in Essence 

Consider the (non-unitary) Schrodinger picture evolution equation 

where H is the usual Hamiltonian, w(t') is an arbitrary function of white noise class, A is a 
Hermitian operator (^|a„) = a„|a„)), A is a collapse rate parameter, T is the time-ordering 
operator and h ^ 1. Associated with this, the probability rule 

t/dt 

Pt{w)Dw =^?A, t\i;, t)u, n dw{t,)/{2TT\/dty/^ (13) 

is defined, which gives the probability that nature chooses a noise which lies in the range 
{w{t'), w{t') + dw{t')} (for calculational purposes, time is discretized, with to = 0). 

Equations (|12|l and l|13|) contain the essential features of CSL, and are all that is needed to 
discuss the simplest collapse behavior. Set H = 0, so there is no competition between collapse 
and the usual Schrodinger evolution, and let the initial statevcctor be 0) — X^nQ^nl^n)- 
Equations H12|) and Ijl^fl become 

\^,t)^ = ^a„|a„)e-(^^)"'-fo'^*'[-(*')-2^''"l^ (14a) 

n 

P,(W) = ^|^j2g-(2A)-VJdt'[»(0-2Aa„]^_ (^4^) 
n 

When the un-normahzed state vector in (|14a|l is divided by (w) and so normahzed, the 
squared amphtudes are 

which are readily shown to satisfy Hlla() . Ijllb|l . and (|llc() in the form x}/'^ {oo)xm'^ (oo) — 
(m ^ n) (which does not change the argument in section 3.1, but makes for an easier calculation). 
Thus, H14a|l . (|14b|) describe collapse dynamics. 
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To describe collapse to a joint eigenstate of a set of mutually commuting operators , replace 
(4A)-i[w(t') - 2Ai]2 in the exponent of (O by Er(4A)"Mw'''(^') " 2Ai'']2. The interaction 
picture state vector in this case is (I12II multiplied by exp(iHt): 

= re-(4^)"/o'it'i:.[-'-(t')-2AA-(0]^l^^o), (15) 

where A''{t') = exp{iHt')A^ exp{-iHt'). The density matrix follows from (^SJ, ifT^ : 

m = J Pt{w)Dw\i>,t)^^{ij,t\/Ptiw) = re-^/2/*dt'i:,[A2(t')-A5,(t'r^(0) (16) 

where (A^(t')) appears to the left (right) of p{0), and is time-ordered (time reverse- 

ordered). In the example described by pil) . the density matrix l|l()|) is 

m = ^e-(^*/2)(«"--")'a„a;,|a„)(a„|, 
which encapsulates the ensemble's collapse behavior. 



3.3 CSL 

The CSL proposal (Pearle 1989) is that collapse is engendered by distinctions between states at 
each point of space, so the index r of in (|15|l becomes x, 

\ij,t)^: = Te-^^^'>'' ^0 1 'it''i-'l-'(^''t')-^^M-',t'}f\^^0)^ (17) 

and the distinction looked at is mass density. However, one cannot make the choice ^(x, 0) = 
M(x), where M(x) = "^iC7 W^i(x) is the mass density operator (m^ is the mass of the ith 
type of particle, so mg, nip, rrin... are the masses of electrons, protons, neutrons..., and if|(x) 
is the creation operator for such a particle at location x), because this entails an infinite rate 
of energy increase of particles f (|23|) with a = 0). Instead, adapting a "gaussian smearing" idea 
from the Ghirardi et al. (1986) spontaneous localization (SL) model (see section 3.6), choose 
A" as, essentially, proportional to the mass in a sphere of radius a about x: 

i(x, t) ^ I dz^e-(2''')"^(--^)%-^^* (18) 

[■Ka^yi'^ J rup 

The parameter value choices of SL, A « lO^^^sec"^ (according to (17), the collapse rate for 
protons) and a w 10~^cm are, so far, consistent with experiment (see section r3.4f) . and will be 
adopted here. 

The density matrix associated with p7|l is, as in (|16|l . 

m ^ Te-(^/2) /* dt'dx'[i,(x',tVl«(x',O]=^(0), (19) 
which satisfies the differential equation 

^^-^J dx'[i(x',i),[i(x',t),p(t)]] (20) 
of Lindblad-Kossakowski form. 



3.4 Consequences of CSL 

Since the state vector dynamics of CSL is different from that of standard quantum theory, there 
are phenomena for which the two make different predictions, allowing for experimental tests. 
Consider an iV- particle system with position operators Xi (X^jx) = x^jx)). Substitution of 
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yi(x') from p8|l in the Schrodinger picture version of (|20|) . integration over x', and utilization 



of 



N 



/(z)Af(z)lx) = ^ m,/(X,)<5(z - XO|x) 





, TV N 

2 ^ TOp mp 




which is a useful form for calculations first suggested by Pearle and Squires in 1994. 
3.4.1 Interference 

Consider the collapse rate of an initial state {(p) — ai|l) + q;2|2), where |1), |2) describe a 
clump of matter, of size << a, at diff'erent locations with separation >> a. Electrons may be 
neglected because of their small collapse rate compared to the much more massive nucleons, and 
the nucleon mass difference may be neglected. In using H21II to calculate d{l\p{t)\2) /dt, since 
exp[— (4a^)^^(Xi ~ X^)^] « 1 when acting on state |1) or |2), and « when acts on |1) and 
Xj acts on |2), (|21|l yields, for N nucleons, the collapse rate XN^: 



If the clump undergoes a two-slit interference experiment, where the size and separation condi- 
tions above are satisfied for ATsec, and if the result agrees with the standard quantum theory 
prediction to 1%, it also agrees with CSL provided A"-'^ > lOOA^^AT. So far, interference exper- 
iments with TV as large as ~ 10^ have been performed by Nairz, Arndt and Zeilinger in 2000. 
The SL value of A would be testable, that is, the quantum predicted interference pattern would 
be "washed out" to 1% accuracy, if the clump were an « 10~^cm radius sphere of mercury, 
which contains N « 10^ nucleons, interfered for AT — .Olsec. Currently envisioned but not 
yet performed experiments (e.g., by Marshall, Simon, Penrose and Bouwmester in 2003) have 
been analyzed (e.g., by Bassi, Ippoliti and Adler in 2004 and by Adler in 2005) involving a su- 
perposition of a larger clump of matter in slightly displaced positions, entangled with a photon 
whose interference pattern is measured: these proposed experiments are still too crude to detect 
the SL value of A, or the gravitationally-based collapse rate proposed by Penrose in 1996 (see 
section 0] and papers by Christian in 1999 and 2005). 

3.4.2 Bound State Excitation 

Collapse narrows wave packets, thereby imparting energy to particles, li H = J2iLi Pi/^nii + 
t^(xi, ...xat), it is straightforward to calculate from H21|l that 



For a nucleon, the mean rate of energy increase is quite small, sa 3 x lO^^^eV/sec. However, 
deviations from the mean can be significantly greater. 

For, (|21() predicts excitation of atoms and nucleii. Let \Eq) be an initial bound energy 
eigenstate. Expanding H21|l in a power series in (bound state size/a)^, the excitation rate of 
state \E) is 



dimm) 



t{mt),H]\2)^xN^{i\p{t)m. 



(22) 



dt 




(23) 



d{E\m\E) 



A 



N N 



T = 



dt 



2a2 



{E\ V !!^X,\Eo)-{Eo\y] ^%\E) + 0(size/a)4. (24) 
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Since \Eq), \E) are eigenstates of the center of mass operator J^iLi X^i^i with eigen- 

value 0, the dipole contribution exphcitly given in l|24() vanishes identicaUy. This leaves the 
quadrupole contribution as the leading term, which is too small to be measured at present. 

However, the choice of A(x) as mass density operator was made only after experimental 
indication. Let gi replace nii/mp in (|21|l . H24|l . so that Xgf is the collapse rate for the «th 
particle. Then, experiments looking for the radiation expected from "spontaneously" excited 
atoms and nucleii, in large amounts of matter for a long time, as shown by CoUett, Pearle, 
Avignone and Nussinov in 1995, Pearle, Ring, Collar and Avignone in 1999 and Jones, Pearle 
and Ring in 2004, have placed the following limits: 

\ge/gp - rrie/mpl < Unie/mp, \gn/gp - mn/mp\ < 3(m„ - mp)/mp. 



3.4.3 Random Walk 

According to (|17(l . (jl^fl . the center of mass wave packet, of a piece of matter of size « a or 
smaller, containing N nucleons, achieves equilibrium size s in a characteristic time Tj, and 
undergoes a random walk through a root mean square distance AQ: 

a^h 1'/* Nnips^ ftAi/2i3/2 

, » AQ» . (25) 



XmpN^ 

The results in H25|l were obtained by CoUett and Pearle in 2003. These quantitative results can 
be qualitatively understood as follows. 

In time At, the usual Schrodinger equation expands a wave packet of size s to w s + 
{h/Nmps)At. CSL collapse, by itself, narrows the wave packet to « s[l — XN'^ {s / a)'^ At] . The 
condition of no change in s is the result quoted above. Ts is the time it takes the Schrodinger 
evolution to expand a wavepacket near size s to size s: {h/NmpS)Ts ~ s. The t^^'^ dependence 
of AQ arises because this is a random walk without damping (unlike Brownian motion, where 
AQ ~ t^/^). The mean energy increase ~ XNfi^m~^a~^t of H23II implies the root-mean-square 
velocity increase ~ [A^j^rn" ^a~^i]^/^, whose product with t is AQ. 

For example, a sphere of density Igm/cc and radius lO^^cm has s « 4 x 10~^cm, Ts ~ 0.6sec 
and AQ « b[t in daysj'^/^cm. At the reported achieved low pressure of 5 x 10~^^Torr at 4.2°K 
reported by Gabrielse's group in 1990, the mean collision time with gas molecules is sa 80min, 
over which AQ « 0.7mm. Thus, observation of this effect should be feasible. 



3.5 Further Remarks 

It is possible to define energy for the io(x, t) field so that total energy is conserved: as the 
particles gain energy, the w-field loses energy, as shown by Pearle in 2005. 

Attempts to construct a special relativistic CSL-type model have not yet succeeded although 
Pearle in 1990, 1992, 1999, Ghirardi, Grassi and Pearle in 1990 and Nicrosini and Rimini in 
2003 have made valiant attempts. The problem is that the white noise field ^(x, t) contains all 
wavelengths and frequencies, exciting the vacuum in lowest order in A to produce particles at 
the unacceptable rate of infinite energy/sec-cc. Collapse models which utilize a colored noise 
field w have a similar problem in higher order. In 2005, Pearle suggested a "quasi-relativistic" 
model which reduces to CSL in the low speed limit. 

CSL is a phenomenological model which describes dynamical collapse so as to achieve S. 
Besides needing decisive experimental verification, it needs identification of the i/7(x, t) field 
with a physical entity. 

Other collapse models which have been investigated are briefiy described below. 



3.6 Spontaneous Localization Model (SL) 

The SL model of Ghirardi et al. (1986), although superseded by CSL, is historically important 
and conceptually valuable. Let H = for simplicity, and consider a single particle whose wave 
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function at time t is ip{'x.,t). Over the next interval dt, with probabihty 1 — \dt, it does not 
change. With probabihty \dt it does change, by being "spontaneously localized" or "hit." A hit 
means that the new (unnormalized) wavefunction suddenly becomes 

^/'(x,f + = V'(x,i)(7ra2)~3/^e-(2'j')"'(x-z)^ with probability Xdtdz j dxlV^x, i + 

Thus z, the "center" of the hit, is most likely to be located where the wavefunction is large. For 
a single particle in the superposition described in section 15.4. II a single hit is overwhelmingly 
likely to reduce the wave function to one or the other location, with total probability |aip, at 
the rate A. 

For an N particle clump, it is considered that each particle has the same independent prob- 
ability, Xdt, of being hit. But, for the example in section [3.4.11 a single hit on any particle in 
one location of the clump has the effect of multiplying the wave function part describing the 
clump in the other location by the tail of the gaussian, thereby collapsing the wave function at 
the rate XN . 

By use of the gaussian hit rather than a delta function hit, SL solves the problem of giving 
too much energy to particles as mentioned in section l53l It also solves the problem of achieving a 
slow collapse rate for a superposition of small objects and a fast collapse rate for a superposition 
of large objects. However, the SL hits on individual particles destroys the (anti-) symmetry 
of wave functions. The CSL collapse toward mass density eigenstates removes that problem. 
Also, while SL modifies the Schrodinger evolution of a wave function, it involves discontinuous 
dynamics and so is not described by a modified Schrodinger equation as is CSL. 



4 Other Models 

For a single (low-energy) particle, the polar decomposition \1/ — i?e^'/''-'"^ of the Schrodinger 
equation implies two real equations, 

(the continuity equation for — \^^) and 

where Q = —{h^/2m)V^R/R is the "quantum potential'. (These equations have an obvious 
generalisation to higher-dimensional configuration space.) In 1926, Madelung proposed that one 
should start from 1)26(1 and l(27|l - regarded as hydrodynamical equations for a classical charged 
fluid with mass density mR^ and fluid velocity VS/m - and construct '5 = ReS^^^^^ from the 
solutions. 

This "hydrodynamical" interpretation suffers from many difficulties, especially for many- 
body systems. In any case, a criticism by Wallstrom (1994) seems decisive: (|26|) and (|27ll (and 
their higher-dimensional analogs) are not, in fact, equivalent to the Schrodinger equation. For, 
as usually understood, the quantum wave function ^' is a single- valued and continuous complex 
field, which typically possesses nodes (^' = 0), in the neighborhood of which the phase S is 
multivalued, with values differing by integral multiples of 27r/i. If one allows S in l|26|l . (|27|l to 
be multi- valued, there is no reason why the allowed values should differ by integral multiples of 
27rfi, and in general ^' will not be single- valued. On the other hand, if one restricts S in ((26|l . 
H27|l to be single- valued, one will exclude wave functions - such as those of nonzero angular 
momentum - with a multivalued phase. (This problem does not exist in pilot-wave theory as 
we have presented it here, where 5* is regarded as a basic entity.) 

Stochastic mechanics, introduced by Fenyes in 1952 and Nelson (1966), has particle trajecto- 
ries x(t) obeying a "forward" stochastic differential equation dx{t) — h{x{t),t)dt + dw{t), where 
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b is a drift (equal to the mean forward velocity) and w a Wiener process, and also a similar 
"backward" equation. Defining the "current velocity" v — i(b + b,), where b* is the mean 
backward velocity, and using an appropriate time-symmetric definition of mean acceleration, 
one may impose a stochastic version of Newton's second law. If one assumes, in addition, that 
V is a gradient (v — VS/m for some S*), then one obtains (|26(l . (|27|l with R = y/p, where p 
is the particle density. Defining 5* = y^e'-'/''^'^, it appears that one recovers the Schrodinger 
equation for the derived quantity 5". However, again, there is no reason why S should have the 
specific multivalued structure required for the phase of a single-valued complex field. It then 
seems that, despite appearances, quantum theory cannot in fact be recovered from stochastic 
mechanics (Wallstrom 1994). The same problem occurs in models that use stochastic mechanics 
as an intermediate step (e.g., Markopoulou and Smolin in 2004): the Schrodinger equation is 
obtained only for exceptional, nodeless wave functions. 

Bohm and Bub (1966) first proposed dynamical wave function collapse through deterministic 
evolution. Their collapse outcome is determined by the value of a Wiener-Siegel hidden variable 
(a variable distributed uniformly over the unit hypersphere in a Hilbert space identical to that 
of the statevector). In 1976, Pearle proposed dynamical wave fmiction collapse equations where 
the collapse outcome is determined by a random variable, and suggested (Pearle 1979) that 
the modified Schrodinger equation be formulated as an Ito stochastic differential equation, a 
suggestion which has been widely followed. (The equation for the state vector given here, which 
is physically more transparent, has its time derivative equivalent to a Stratonovich stochastic 
differential equation, which is readily converted to the Ito form.) The importance of having 
the density matrix describing collapse be of the Lindblad-Kossakowski form was emphasized by 
Gisin in 1984 and Diosi in 1988. The stochastic differential Schrodinger equation which achieves 
this was found independently by Diosi in 1988 and by Belavkin, Gisin and Pearle in separate 
papers in 1989 (see also Ghirardi et al. 1990). 

A gravitationally motivated stochastic collapse dynamics was proposed by Diosi in 1989 (and 
somewhat corrected by Ghirardi et al. in 1990). Penrose emphasized in 1996 that a quantum 
state, such as that describing a mass in a superposition of two places, puts the associated space- 
time geometry also in a superposition, and has argued that this should lead to wave function 
collapse. He suggests that the collapse time should be ^ /i/Ai?, where Ai? is the gravitational 
potential energy change obtained by actually displacing two such masses: for example, the 
collapse time « H/{Gm?/R), where the mass is to, its size is R, and the displacement is ~ i? or 
larger. No specific dynamics is offered, just the vision that this will be a property of a correct 
future quantum theory of gravity. 

Collapse to energy eigenstates was first proposed by Bedford and Wang in 1975 and 1977 
and, in the context of stochastic collapse (e.g., (11) with A = H), by Milburn in 1991 and 
Hughston in 1996, but it has been argued by Finkelstein in 1993 and Pearle in 2004 that such 
energy-driven collapse cannot give a satisfactory picture of the macroscopic world. Percival in 
1995 and in a 1998 book, and Fivel in 1997 have discussed energy-driven collapse for microscopic 
situations. 

Adler (2004) has presented a classical theory (a hidden-variables theory) from which it is 
argued that quantum theory "emerges" at the ensemble level. The classical variables are N x N 
matrix field amplitudes at points of space. They obey appropriate classical Hamiltonian dy- 
namical equations which he calls "trace dynamics," since the Hamiltonian, Lagrangian, Poisson 
bracket, etc. expressions have the form of the trace of products of matrices and their sums with 
constant coefficients. Using classical statistical mechanics, canonical ensemble averages of (suit- 
ably projected) products of fields are analyzed and it is argued that they obey all the properties 
associated with Wightman functions, from which quantum field theory, and its non-relativistic 
limit quantum mechanics, may be derived. As well as obtaining the algebra of quantum theory 
in this way, it is argued that statistical fiuctuations around the canonical ensemble can give rise 
to wave function collapse behavior, of the kind discussed here, both energy-driven and CSL-type 
mass-density-driven collapse. The Hamiltonian needed for this theory to work is not provided 
but, as the argument progresses, its necessary features are delimited. 
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